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LARS Contract Report 110984 


Introd u ction 

This report covers work carried out on Landsat-4/5 data 
quality analysis under NASA Contract NAS5-26859 for the per- 
iod August 10 throu^ November 9» 1984. Focus was on esti- 
mation of two-dimensional point-spread function estimation. 

A brief description is included in the results section. 
Preparation began on the final report for the project. 

ProblCTis 

Landsat-5 data needed for the final phase of the study 
were not received during the quarter, therefore no analyses 
could be done on evaluation of these data. Landsat-5 TM 
data were received for three sites on November 15, 1984. 
There is approximately one quarter of funding left in the 
project, so all work will be done within these remaining 
resources . 


Publications 

No publications were generated during the quarter. 
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Recommendations 

No recommendations are made In this report. 


Funds Expended 

The funds expended In the project are reported periodi- 
cally by the Purdue University Office of Contract and Grant 
Business Affairs to the sponsor on NASA Fonn 533M. ^eclflc 
disclosure of funds expended In this report Is not a policy 
of the University. 


Significant Results 

One of the system parameters estimated under this 
NASA-sponsored program has been the PSF or blurring function 
of the Landsat Thematic Mapper. The PSF Is the result of 
the blurring effect of the atmosphere, electronics, optics, 
and data-processlng operations that take place on the sig- 
nal. With a good estimate of the blurring function on an 
imaging system, the next logical step Is to use this Infor- 
mation to restore the data radlometrlcally. 

A slgnficant amount of work has been done on the area 
of image restoration. The classical approach that has 
mostly been followed is in the direction of the Wiener fil- 
ter theory. There are major reasons against the usage of 
Wiener filters on Landsat data: 
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1. It la assvmed that the image is a stationary 
random field » which is not a good assianption 
for agricultural fields. 

2. The filter design requires statistics of the 
input signal. 

3. A (fiiSE criterion is usedy which has been shown 
not to be optimal for Images whose final pur- 
pose is to be analyzed by a human observer. 

Because of these limitations of the Wiener filter » a 
different approach was used which is base! on minimizing a 
narrowness measure in the composite PSF (CPSF*)| subject 
to constraints on the noise levels at the output of the fil- 
ter (ChUt Stuller). Some of the advantages of this method 
are: 

1. The resulting discrete filter is signal inde- 
pendent. The only a-priori information needed 
in the design is the PSF and the autocorrela- 
tion function of the noise. 

2. The restoration filter is a FIR filter whose 
duration can be conveniently selected (within 
a certain range). 

• The composite point-spread function (CPSF) is defined as 
the convolution of PSF and the restoration filter. 



We have been looking at new methods of introducing full 
control over the sidelobe levels present on the CPSF. The 
drawback of introducing constraints on the sidelobe levels 
is that an analytical solution for the filter coefficients 
cannot be obtained. But, by properly formulating the prob- 
lem, this can easily be solved using nonlinear programming 
techniques. In particular, we are using a generalized 
reduced-gradient algorit'nm. For the case of space invariant 
PSF and stationary noise, the filter design has to be per- 
formed only once, making it very attractive for processing 
large volumes of data, as in the case of Landsat. 

We are presently evaluating the performance of this 
image-restoration technique on Thematic Mapper data from the 
Iowa and Chicago sites. Also, we plan to combine th ^ 
restoration filter with an interpolating function to simul- 
taneously obtain improved radiometric and spatial resolu- 
tion. 
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